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Lattice QCD can contribute to the search for new physics in b — > s 
decays by providing first-principle calculations of B — > K^ form fac- 
tors. Preliminary results are presented here which complement sum rule 
determinations by being done at large q 2 and which improve upon pre- 
vious lattice calculations by working directly in the physical b sector on 
unquenched gauge field configurations. 
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1 Introduction 



The b — > s flavour changing neutral current transition is suppressed in the standard 
model. Dominant contributions to rare B decays such as B — >■ K^l + l~ come from 
loop diagrams: box and penguin diagrams. These rare B decays are good windows 
for looking for new physics: New particles beyond the standard model could appear 
in the loops and change the decay widths of those rare decays. 

The starting point of theoretical calculations of weak decays of hadrons is the 
effective weak Hamiltonian. In the standard model, there are ten operators in the 
effective Hamiltonian for radiative and semileptonic decays. The dominant short 
distance contributions are from effective local operators Q 7 , Q 9 and Q w , which come 
from the penguin and box diagrams. 

In quantum chromo dynamics (QCD), quarks are confined in color singlets. The 
b — > s transition happens inside hadrons. Therefore the matrix elements of the 
above three local operators have to be computed using non-perturbative methods, for 
example, lattice QCD. Those matrix elements can be parametrized by form factors 
according to their Lorentz structures. In total, there are ten form factors for the 
quark currents in Q 7 , Q 9 , and Q w , and our aim is to calculate these on the lattice 
with dynamical simulations. 

More details of our calculation strategy and our definitions of the form factors 
can be found in Ref. [1]. Here we update our progress in the extraction of the form 
factors. In the next section, we present our lattice setup and then we show some 
preliminary results in the last section. 

2 Lattice setup 

We use configurations from the MIMhQ Lattice Computation (MILC) Collaboration, 
which are 2 + 1 flavour dynamical simulations using 0(a 2 ) and tadpole improved 
staggered fermions (AsqTad) [2]. We currently have data from two lattice spacings. 



a(fm) 




Volume 


N t x N 

1 > con f ™ 1 v src 


am val 


coarse ~0.12 


0.007/0.05 
0.02/0.05 


20 3 x 64 
20 3 x 64 


2109 x 8 
2052 x 8 


0.007/0.04 
0.02/0.04 


fine ~0.09 


0.0062/0.031 


28 3 x 96 


1910 x 8 


0.0062/0.031 



Table 1: Parameters of lattices being used in this study. N src is the number of point 
sources used on each configuration. 

At the coarse lattice spacing we have two different light quark masses. The lightest 
* Multiple Instruction stream, Multiple Data stream. 
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quark mass gives a pion mass of about 300 MeV. The parameters of our calculation are 
collected in Table [TJ On each configuration, eight point sources are used to increase 
statistics. In [lj, we found Z2 x Z2 random wall sources were inefficient for vector 
mesons and heavy-light mesons in reducing statistical errors of correlators (random 
wall source methods allow one to approximately obtain all to all correlators and thus 
possibly to reduce statistical errors). Therefore we now use several point sources. 

The light valance quarks are also AsqTad fermions as the sea quarks. For the 
heavy b quark, we use the (moving-)non-relativistic QCD (NRQCD) action [3], which 
is expanded up to and including 0(AQ CD /ml). We can work directly at the physical 
b quark mass, thus no extrapolation up to is required. In the lattice heavy-light 
currents, the heavy quark expansion includes order l/m&. 

We compute 2-point functions for the heavy-light B meson and the light-light 
final state mesons as well as 3-point functions with the current operators inserted. 
Then we fit these correlation functions with the Bayesian fitting method [3]. From 
the fitted ground state energies and amplitudes, one can extract the matrix elements 
of the current operators and then the form factors. The detailed formulas can be 
found in Refs. [H [5]. 

Matching factors of the lattice current operators to the MS scheme were computed 
perturbative in Ref. (6j to one loop. We set a s = 0.3 below to get the values of these 
matching factors. 

Previous lattice calculations of the above form factors are all quenched calcula- 
tions. And an extrapolation up to the physical b quark mass is needed for the heavy 
quark. See, e.g., [7j and the references therein. 

3 Preliminary results 

We tried both separate fitting and simultaneous fitting of the 2- and 3-point functions. 
The simultaneous fitting, where common energy parameters are used, gives better 
results with smaller statistical errors. Examples are shown in Fig. [1] for the ground 
state amplitude in a 3-point function for B — > K and the B meson energy (in lattice 
units) obtained from different fits. Note a mass shifting term needs to be added for 
the B meson energy due to the use of (moving-)NRQCD. In Fig. [TJthe horizontal axis 
is the number of (ground and excited) states in the Bayesian fitting functions. Equal 
number of opposite parity states are also in the fitting functions due to the use of 
staggered fermions. The fitting results stabilize when more than 6 normal states and 
6 opposite parity states are used. 

Preliminary results of some form factors for B — > K and B — > K* are given in 
Fig. [2] and Fig. [3] respectively. They are from the coarse lattice with light valance 
quark masses 0.007/0.04. The l/m& corrections to the currents are not included yet, 
but will be soon. A first fitting to the 3-point functions shows that these corrections 
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Figure 1: Comparison of separate Bayesian fitting and simultaneous Bayesian fitting 
to heavy-light 2-point (hl2pt) and heavy-light 3-point (hl3pt) functions for the ground 
state amplitude in a 3-point function (left graph) and the B meson energy aE 1 ^ 1 (right 
graph) in lattice units. The horizontal axis is the number of (ground and excited) 
states in the fitting functions. 
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Figure 2: Preliminary results for the 
form factors fo, f + and fx for B — > 
K decays, obtained from simultaneous 
Bayesian fits. The left-most points have 
v = (0,0,0), p B = and p K = 2n/L ■ 
(-2,0,0). 



Figure 3: Preliminary results for the form 
factor Ti and T 2 for B — > K* decays, ob- 
tained from simultaneous Bayesian fits. 
The left-most points have v = (0,0,0), 
p B = and p K . = 2ix / L ■ (-2, 0, 0). 
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are small. For the other quark mass and for the fine lattice spacing, data analysis is 
going on. We will have more data points at lower q 2 after we analyze the correlators 
with non-zero velocity v for the B meson in the moving NRQCD action. 

For the radiative decay B — > jK*, we want the form factor T(0) = Ti(0) = ^2(0) 
to get the branching fraction. To extrapolate our results to the q 2 = limit, we need 
some ansatz. In Fig. H] we do an extrapolation using the B-K ansatz [8] 
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Figure 4: Extrapolation of T\ and T 2 to q 2 = using the B-K ansatz in Eq.flT]). 

The statistical errors will be reduced when more data points at low q 2 are included. 

Eventually we need to extrapolate to the physical pion mass and the continuum 
limit to really compare with experiments. In a recent paper by H. Na et al. [9] on 
D — >■ K form factors, the extrapolation in q 2 , light quark mass and lattice spacing 
were done all together using a "simultaneous modified ^-expansion" . We may try this 
method in the future. 
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